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Abstract 

In earlier Letters, we adopted a complex approach to quantum processes in the formation 
and evaporation of black holes. Taking Feynman's + ie prescription, rather than one of 
the more usual approaches, we calculated the quantum amplitude (not just the probability 
density) for final weak-field configurations following gravitational collapse to a black hole 
with subsequent evaporation. What we have done is to find quantum amplitudes relating 
to a pure state at late times following black-hole matter collapse. Such pure states are 
then shown to be susceptible to a description in terms of coherent and squeezed states - in 
practice, this description is not very different from that for the well-known highly-squeezed 
final state of the relic radiation background in inflationary cosmology. The simplest such 
collapse model involves Einstein gravity with a massless scalar field. The Feynman ap- 
proach involves making the boundary-value problem for gravity and a massless scalar 
field well-posed. To define this, let T be the proper-time separation, measured at spatial 
infinity, between two space-like hypersurfaces on which initial (collapse) and final (evap- 
oration) data are posed. Then, in this approach, one rotates T — > |T|exp(— id) into the 
lower half-plane. In an adiabatic approximation, the resulting quantum amplitude may be 
expressed in terms of generalised coherent states of the quantum oscillator, and a physical 
interpretation is given. A squeezed-state representation, as above, then follows. 

1. Introduction 

We begin by describing Feynman's + ie approach [1] in the context of black-hole 
evaporation. In [2-12], this treatment was described and applied to the calculation of 
quantum amplitudes (not just probabilities) for particle production, following gravitational 
collapse to a black hole. Suppose, for definiteness, that one's Lagrangian contains Einstein 
gravity coupled to a real massless scalar field. Asymptotically-flat initial data are posed 
on an initial space-like hypersurface E/ , and final data on a surface E p , separated from 
Tij by a (large) real Lorentzian time-interval T , as measured at spatial infinity. Suppose 
further, for simplicity, that the initial data on E/ are spherically symmetric, corresponding 
to a diffuse slowly-moving initial matter distribution. The final data for gravity + scalar 
are taken to have a 'background' spherically-symmetric part, plus small non-spherical 
perturbations, which correspond to gravitons and massless-scalar particles. 

Following Feynman's +ie procedure [1], one rotates the time-interval T into the com- 
plex: T — > \T\ exp(— id) , with < 5 < n/2 . The classical boundary- value problem, for a 
complex 4-metric g^ v and scalar field <p given the above data on E/ , Tip , is then expected 
to be well-posed, unlike the ill-posed case 6 = (or equivalently T real) [3,13,14]. One 

(2) 

can evaluate the second-variation classical action S^iass as a functional of the (still real) 
boundary data and as a function of the complex variable T . One then computes the cor- 
responding semi-classical quantum amplitude, proportional to exp^S^j^J , and can also 
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include loop corrections, if appropriate. Finally, the Lorentzian quantum amplitude for 
black-hole evaporation (again, not just the probability density) is recovered by taking the 
limit as 5 — > + . 

In this Letter, we study such black-hole evaporation amplitudes, which were con- 
structed in detail in [5-8, 10-12], but now in the context of coherent states [15], which re- 
semble 'classical states', and of squeezed states [16], which are purely quantum-mechanical. 
Although our motivation originated with the question of black-hole radiation, there are 
also strong connections between this work and the study of the relic Cosmic Microwave 
Background Radiation (CMBR) induced by inflationary cosmological perturbations. 

In inflationary cosmology, the field modes are in their adiabatic ground state, with 
short wavelengths near the start of inflation. Due to the accelerated expansion of the 
universe during inflation, quantum fluctuations are amplified into macroscopic or classical 
perturbations. The early-time fluctuations lead to the formation of large-scale structure 
in the universe, and also contribute to the anisotropies in the CMBR. The final state 
for the perturbations is a two-mode highly-squeezed state for modes whose radius is much 
greater than the Hubble radius [17], pairs of field quanta being produced at late times with 
opposite momenta. Tensor (s = 2) fluctuations in the metric, for example, are predicted 
to give rise to relic gravitational waves. By comparison, electromagnetic waves (s = 1) 
cannot be squeezed in the same way. 

In either case, cosmological or black-hole, one works within an adiabatic approxi- 
mation for the perturbative modes. Writing k for a typical perturbative frequency, one 
requires k^>H in the cosmological case, where H = (a/a) and a(t) is the scale factor. 
In the black-hole case, the space-time geometry at late times, in the region containing a 
stream of outgoing radiation, is given by a Vaidya metric [4,8,18,19] with a slowly- varying 
'mass function' m(t, r) . The adiabatic condition then reads k 3> \m/m\ . 

In applying the squeezed-state formalism, one finds, in the case of cosmological pertur- 
bations, that these evolve essentially according to a set of Schrodinger equations [20]. Such 
perturbations, whether of density, rotational or gravitational type, starting in an initial 
vacuum state, are transformed into a highly-squeezed vacuum state, with many particles, 
having a large variance in their amplitude (particle number), but small (squeezed) phase 
variations. The squeezing of cosmological perturbations may be suppressed at small wave- 
lengths, but it should be present at long wavelengths, especially for gravitational waves 
[21]. These perturbations also induce the anisotropies at large angular scales, as observed 
in the CMBR. Their wavelengths today are comparable with or greater than the Hubble 
radius. The above amplification of the initial zero-point fluctuations gives rise to stand- 
ing waves with a fixed phase, rather than travelling waves. The relic perturbations in 
the high-squeezing or WKB limit can be described as a stochastic collection of stand- 
ing waves. Although this paragraph has reviewed the application to cosmology, a similar 
picture emerges in the application to black-hole evaporation. 

Sec. 2 outlines the main features of the above complex approach to the calculation 
of quantum amplitudes (not just probabilities) for data (spins s = , 1 , 2) prescribed on 
a late-time final hypersurface T,p . This requires a rotation: T — > |T| exp(— iS) into the 
lower half-plane. The resulting amplitudes are then related to coherent and squeezed 
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states. Sees. 3, 4, 5 describe coherent states, generalised coherent states and squeezed states, 
respectively. In Sec.6, the small angle 8 (above), through which the time T at infinity 
is rotated into the complex, is related to the large amount of squeezing which has been 
applied to give the final state. Sec. 7 contains a brief Conclusion. 

2. The quantum amplitude for late-time data 

Consider first the case of a rotation into the complex of the time-interval T , measured 
at spatial infinity, by a moderately small angle 5 , as above. One expects that the resulting 
classical solution (g^ , 0) of the coupled Einstein/massless-scalar field equations is slightly 
complexified, by comparison with a Lorentzian-signature solution. By suitable choice of 
coordinates (t,r,0,<p), the spherically-symmetric 'background' part of the metric may be 
written in the form [2,6] 

ds 2 = - e b dt 2 + e a dr 2 + r 2 (d9 2 + sm 2 9d^ 2 ) , (2.1) 

where b = b(t , r) , a = a(t , r) , and the spherically-symmetric 'background' part $ of the 
scalar field has the form $ = $(£,r) . The coupled Lorentzian-signature Einstein/scalar 
field equations for this spherically-symmetric configuration are given by the analytic con- 
tinuation of the Riemannian field equations Eqs. (3.7-11) of [5], on making the replacement 

t = rexp(-itf) , (2.2) 

where r is the 'Riemannian time-coordinate' of [5], and where the real number should 
be rotated from to n/2 . 

Small non-spherical perturbations in the boundary data given on the final late-time 
hypersurface E p consist of the perturbed part of the intrinsic 3-dimensional spatial metric 
hijF on Ef, together with the perturbations in the scalar field on Sj? . As above, these 
correspond to gravitons and to massless-scalar particles, propagating on the spherically- 
symmetric classical background [g^ , <E>) . For example, the linearised scalar perturbations 
</>W , given [2] by = $ + </>W , may be first decomposed as in Eq.(6) of [2], namely as: 

oo m=£ 

<jP-\t,r,0,tp) = -J2 E Yem(ty Remit, r) . (2.3) 

£=0 m=-£ 

Here, Yi m (Q) denotes the (£,m) spherical harmonic of [22]. The scalar field equation 
decouples for each (£,m), leading to the mode equation 

(>- a )/ 2 drj ''R im - (d t ) 2 R lm - \ (dt(a-b)) (d t R lm j - V £ (t,r)R £rn = , (2.4) 
where 

V,( t ,r) = ^(/(/ + l) + *^) (,5> 
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is real and positive in the Lorentzian-signature case. The 'mass function' m(t , r) , which 
would equal the constant mass M for an exact Schwarzschild geometry [23] , is defined by 



e -(*,r) = 1 _ 2m(*,r) 

r 

An analogous harmonic decomposition can be given for weak gravitational-wave perturba- 
tions about the spherical background [10]. 

In most regions of the classical space-time, except for the central region where the black 
hole is formed, the metric functions a(t,r) and b(t,r) vary 'slowly' or 'adiabatically'. In 
this case, one can consider a radial mode solution for (say) a perturbed scalar field, of the 
form [6] 

Rem(t, r) ~ exp(ikt) £k£m{t,r) , (2.7) 

where £k£m(t , v) varies 'slowly' with respect to t . This will occur near spatial infinity, and 
it will also occur, provided that the time-interval T is sufficiently large, in a neighbourhood 
of the final hypersurface T,p . The mode equation (2.4,5) then reduces [6] to 

e^ 2 I [e^' 2 9 -^j + (k 2 - V t ) ium = . (2.8) 

The spherically-symmetric background metric in this region can be represented to high ac- 
curacy by a Vaidya metric [8,18,19], which describes the (on average) spherically-symmetric 
outflow of massless matter. The principal condition for the validity of the adiabatic ex- 
pansion is [6] that 

\k\ > \m/m\ . (2.9) 

In analysing the behaviour of the radial mode equation (2.8), it is natural to define 
a generalisation r* of the standard Regge- Wheeler or 'tortoise' coordinate r£ for the 
Schwarzschild geometry [23], according to 

° = e^/ 2 I- . (2.10) 



dr* dr 
The approximate (adiabatic) mode equation (2.8) then reads 



d 2 Ck£m 

dr 



*2 



+ 



(k 2 - V^Mm = . (2.11) 



We consider here, for definiteness, a set of suitable radial functions {(HmW} on the 
final surface E p , since it is here that the non-trivial boundary data are posed. Since 
the mode equation (2.11) does not depend on the quantum number m, we may choose 
£k£m( r ) = 6cK r ) 5 independently of m. The boundary condition of regularity at the 
spatial origin {r = 0} [6] implies that 

e+1 + n((hA l ^\ 



£ M {r) = constant x (At) " 1 " + O^(Ar) J (2.12) 



4 



as r — > 0_|_ . For the boundary condition on the Cki( r ) as r — > oo , note that the potential 
Vg(r) decreases sufficiently rapidly, as r — > oo , that a real solution to Eq.(2.11) behaves 
near {r = 00} according to 

f fc *(r) ~ (z ki esxpfikr^j + z* kl exp^-ifer^ . (2.13) 

Here, the ^ are certain dimensionless complex coefficients, which must be determined by 
using the differential equation (2.11) together with the regularity conditions. Further [6], 
there is a natural normalisation of the basis {£ki(r)} of radial wave- functions. 

We continue, for purposes of exposition, to study the case of scalar perturbations, 
with a slightly complexified time-interval at infinity, T = \T\ exp(— id) , for < 5 < 7r/2 . 
The relevant boundary data for anisotropic perturbations (p^ of the scalar field 4>f on 
Tip can be described [6] by expanding out the interior classical boundary- value solution 
near E_p in the form 

= - E E r ^ a kirn^t(t,r)^j^rY lm (^) ■ (2.14) 

r i= G m =-l J -°° sm{Kl) 

Here, the real quantities {a k £ m } characterise the final data. 

More generally, for perturbative boundary data for a field of any spin, posed on Ep 
in describing a final state resulting from black-hole evaporation, we denote by {a s kimp} a 
set of analogous 'Fourier-like' coefficients, where s gives the particle spin, k the frequency, 
(£,m) the angular quantum numbers, and P = ±1 the parity (for s^O). For massless 
perturbations of spins s = 0, 1 ,2 [2,3,5-7,10,11], we found that the quantum amplitude 
or wave functional is of semi-classical form, being given by 



{a sk imp}]T = N exp^Sdass {a skim p}; T ) , (2.15) 



where the pre-factor iV depends only on T . Here, S^ass denotes the (second- variation) 
action of the classical infilling solution, as a functional of the boundary data. For simplicity, 
we denote the collection a sk tmP of indices by j . Further, we write Mj for the total (time- 
independent) ADM (Arnowitt-Deser-Misner) mass of the 'space-time', as measured at 
spatial infinity [23] . The ADM mass Mj , which is the limit at large radius of the variable 
mass m(t , r) of the Vaidya metric, is a functional of the final field configurations {aj} on 
£ f , since it depends on the full gravitational field which results from classical solution of 
the complexified boundary-value problem. 

As was found (for example) in the scalar case s = in [2,3,6], the classical action is 
dominated by contributions from frequencies k with the values 

k = k n = — ; n=l,2,3,... . (2.16) 

We also define Akj to be the spacing between neighbouring fcj-values: 

Akj = £ . (2.17) 
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Following [2,3,5-7,10,11], the classical action functional S c \ ass is found to be a sum 
over individual 'harmonics' labelled by j , which depend on the corresponding indices 
{s kj Im P} through the quantity \Aj\ , defined by 



I A, 



2 (-l) s c s 



(£ + s)\ 



Ml 



dj + (~1) S Pa 8 ,-kjlmP 



(2.18) 



Here, the coefficients c s for bosonic spins s are given by Co = 2n , c\ = \ , c 2 = | . The 
quantities Zj are the complex numbers appearing in Eq.(2.13), which arise in solving the 
adiabatic radial mode equation (2.11). This leads to the form of the quantum amplitude: 



{Aj}; t\ = N e~ tMlT ]J V(Aj ; T) 



(2.19) 



where N also depends only on T . 

Taking the classical action iS c i ass in the form found in [6] for the scalar s = case (for 
example), one deduces that the wave functional for given boundary data can be written: 



{Aj}:T = N 



-iMiT 



n 



2i s'm(kjT) 



exp 



- (Akj) kj {Ajl 2 cot(^T) 



. (2.20) 



This will be related to the coherent-state description in the following Section 3. 

3. Coherent States 

It is possible to rewrite the quantum amplitude (2.20) with the help of the Laguerre 
polynomials [24]. First, we introduce the associated Laguerre polynomials L^ m (x) , 
defined by 



for m > n > . The Laguerre polynomials L n (x) [24] are given by 

L n {x) = L^\x) . 
The set {L n (x)} obeys the completeness relation 

oo 

E e " (x/2) e~ iy/2) L n (y) = 5(x,y) 

n=0 



(3.1) 



(3.2) 



(3.3) 



Writing z = x + iy , consider now the function L n (\z\ 2 ) , which appears in Eq.(3.5) below. 
For n > , this cannot be written as a product of two (decoupled) wave functions of x and 
y in an excited state, due to pair correlations [25]. But, in terms of Hermite polynomials 
H p (x) [24], one can write 



L n (x 2 + y z 



(~l) n 

2 2n n ] 



E 

6 



n 



P 



H 2p (x) H2n-2p (y) ■ 



(3.4) 



From this, one can further decompose the quantum amplitude (2.20) as 



oo 

9 '{A s };T\ = Ne~^ T exp(-E, (Akjkj \A d \ 2 /2) ]J £ e" 2 ^ T L n [jfc,- (AJfe,-) \A 



j n=0 

(3.5) 

where E n = (n + |) kj is the quantum energy of the linear harmonic oscillator. Note also 
the dependence of the quantum amplitude on \Aj\ ■ - it is spherically symmetric. 

The Schrodinger-picture wave functions 

<M*;,T) = ^e-^/ 2 ) e" 2 ^ T L n (x,) (3.6) 

2 

appear in the wave-function (3.5), with Xj = kj (Akj) \Aj\ . The wave functions (3.6) 
have a strong connection with the exact solution of the forced-harmonic-oscillator problem 
[26], with Hamiltonian 

H = ^ + + Q F (t) > ( 3 - 7 ) 

where F(t) denotes an external force, n the oscillator mass and uj the oscillator frequency. 
Assume that F(t) = for t < t and for t > T, so that the asymptotic states, at early 
and late times t , are free-oscillator states. One can calculate the amplitude A^m to make 
a transition from the free-oscillator state \m > (with m particles) at early times t < to ? 
to the free-oscillator state \k > at late times t > T . Define the 'Fourier transform' of the 
force: 

(3=[ dtF(t)e~ lujt , (3.8) 

Jt 

and set 

(3.9) 



2 fj,u 

It has been shown [27-29], in the case m>k, that 

^ = ^-<^g) 4 yy m "V^ , (3.10) 

where A is a real phase. This expression also gives A^ m for m < k , since Ak m = A m k is 
symmetric. 

In the adiabatic limit, in which the force F(t) changes extremely slowly, one has z 1 , 
and from general considerations a state which begins as \k > must end up in the same 
state \k > after the time-dependent force has been removed. From Eq.(3.10), one has 



i-kk 



e lX e" (2/2) L k (z) . (3.11) 



The corresponding probability that there should be no change in the number of particles 
is |Afcjt| 2 = e~ z [Lk(z)] 2 . Apart from the introduction of mode labels j denoting the 
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'quantum numbers' {skimP}, together with a necessary re-interpretation for z , these am- 
plitudes are effectively the wave functions (3.5) derived from our boundary-value problem. 

One further viewpoint can be brought to bear on Eq.(3.10), arising from the coherent- 
state representation. Coherent states \a > can be regarded as displaced vacuum states; 
that is, [15] 

\a> = D(a) |0 > , (3.12) 

where 

D(a) = expect a) — a*aj (3.13) 
is a unitary displacement operator, obeying 

D\a) = D- X {a) = D{- a) , (3.14) 

and where the states \a > are eigenstates of the annihilation operator a with complex 
eigenvalue a . Among quantum states for the harmonic oscillator, they are the closest to 
classical states, in that they attain the minimum demanded by the uncertainty principle. 
Coherent states form an over-complete set, and are not orthogonal. In terms of the Fock- 
number eigenstates 

(a^) n 

\n> = ^jL- |0> , (3.15) 
vn! 



one has [29] 

oo 

\a> = e-l^l 2 / 2 Y ^= \n> . (3.16) 



n 



n=0 

The coherent state labelled by a = is the ground state of the oscillator. If, for example, 
the system started in a vacuum state, the amplitude to find it subsequently in a coherent 
state \a > is 

<0|a> = <0\D(a)\0> = e~ |q|2/2 , (3.17) 

up to a phase. 

To make complete contact with the amplitude (3.10), using coherent-state methods, 
we note that, in terms of the displacement operators D(£) : 



<m\D(£)\a> = -yL (£ + a) m exp 



\(\*\ 2 + \z,\ 2 + n*o) 



(3.18) 



and 



oo „ 

or 



<m\D(0\a> = e -l Q l /2 y ^= <m\D(0\n> . (3.19) 



n=0 

On equating these, one finds that 



oo 

2 



n=0 




(l + y) m e~ y ^ = e l?l /2 ^ J- £ n ~ m y n <m\D{t)\n> . (3.20) 



But, from the generating function for the associated Laguerre polynomials [25], 

oo 

(l + y )m e -vx = J- LiT- n \x) y n , \y\ < 1 , (3.21) 

one deduces that the matrix element between initial and final states is 

i 

<m\D{i)\n> = (^j) 2 r~ n e-l^l 2 / 2 4— )(|C| 2 ) , (3.22) 



which agrees with Eq.(3.10), up to an unimportant phase factor. 

4. Generalised coherent states 

These amplitudes can also be interpreted in terms of generalised coherent states of 
the harmonic oscillator [281. Define: 



| ra , a > 

Then, in the Fock representation, 



-iE n t 



D(a(t)) \n > 



(4.1) 



\n,a> = < m\D(a(0))\n > \m > e 



-iE m t 



(4.2) 



m=0 



For generalised coherent states, the ground state (n = 0) is a coherent state and not a 
vacuum state. Generalised coherent states are to the coherent states what the Fock states 
| n > are to the vacuum state, that is, excited coherent states. In addition, denoting by / 
the identity operator, one finds that [27]: 



1 



TV 



d a \n , a > < n , a\ , 



<n,0\n,a> = L n (\a - (3\ 2 ) exp(/T a - 1 (|a| 2 + |/3| 2 )) 



(4.3) 
(4.4) 



e -|/?| 2 /2 



ipi / 2 r / \ 2 

<n,(3\i)> = d 2 a L n (\a-l3\ 2 J e p ' a e" |a| /2 <n,o#> ,(4.5) 

for an arbitrary state \ip > , with the definition: 



J d 2 a = J d 



Re(a) d Im(a) 



(4.6) 



In particular, from Eq.(4.4) with f3 = , one has 

<n,0|n,a> = <n|n,a> =e~' a '/ 2 L n ^|a| 2 j, 

9 



(4.7) 



again giving Eq.(3.6) up to a phase. The initial state should be seen not as a vacuum 
state, but as a Fock state, while the final state should be seen as a generalised coherent 
state. 

As shown by Hollenhorst [30], the amplitudes of Eqs.(3.22) have yet a further inter- 
pretation: they are the matrix elements for a transition from state \k > to state \m > 
under the influence of a linearised gravitational wave, with the force F(t) proportional to 
the Riemann curvature-tensor component Rtxtxif) '■ 

F{t) = ixlR^tit) = - \^l{d t ) 2 h T x l , (4.8) 

where I is the distance between two particles along the x-axis, each being of mass (/i/2) , 
while ti^T is the transverse-traceless gravitational- wave component of the metric [23] , and 
x is the change in the separation of the masses. 

In the context of black- hole evaporation, one expects that the role of the force is played 
by the time-dependent background space-time - which approximates a Vaidya space-time 
in the high-frequency limit at late times [4,8,18,19]. 

An important point which we should mention is that, under the influence of a time- 
dependent force, an initial vacuum state transforms into a coherent state. Below, we discuss 
how, by changing a phase parameter of the perturbations appearing in their frequencies 
(parametric amplification), an initial vacuum state transforms into a squeezed vacuum 
state. This phase is not an oscillator phase, but a small angle, 6 , through which the time 
T at infinity is rotated into the lower complex plane. 

5. Squeezed-state formalism 

In this Section and in the following Sec. 6, we shall see how, by rotating the asymptotic 
Lorentzian time T into the complex plane, and in the case of spherically-symmetric initial 
matter and gravitational fields, one obtains a quantum-mechanical highly-squeezed-state 
interpretation for the final state in black-hole evaporation, in the limit of an infinitesimal 
rotation angle. 

Grishchuk and Sidorov [17] were the first to formulate particle creation in strong 
gravitational fields explicitly in terms of squeezed states, although the formalism does 
appear in Parker's original paper on cosmological particle production [31]. In [17], it 
was shown that relic gravitons (as well as other perturbations), created from zero-point 
quantum fluctuations as the universe evolves, should now be in a strongly squeezed state. 
Squeezing is just the quantum process corresponding to parametric amplification. 

Black-hole radiation in the squeezed-state representation was first discussed in [17]. 
The 'squeeze parameter' rj (see below) was there related to the frequency ujj and the 
black-hole mass M through 

tanh(r.,) = exp(-47rMu;^ . (5.1) 

In this language, the vacuum quantum state in a black-hole space-time for each mode is 
a two-mode squeezed vacuum. However, our approach to squeezed states in black-hole 
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evaporation is new; arising from a two-surface boundary-value problem and Feynman's 
+ ie prescription [1] . We now give a brief account of quantum- mechanical squeezed states. 

A general one-mode squeezed state (or squeezed coherent state) is defined [16] as 



a,z > = £>( 7 ) S(r,<f>) |0 > 



£>(7) S(z) |0 > 



(5.2) 



Here, D{^) is the single-mode displacement operator, and 




)) 



(5.3) 



in terms of annihilation and creation operators a and a< , respectively, together with the 
relation 




SHz)S(z) 



S(z)SHz) 



= 1 



(5.5) 



with 7 given by 



7 



a cosh r + a* e 



2l4> sinh r 



(5.6) 



The state Eq.(5.2) is a Gaussian wave-packet, displaced from the origin in position and 
momentum space. While the (real) squeezing parameter r (0 < r < oo) determines the 
magnitude of the squeezing, the squeezing angle (p (1^1 < 7r /2) gives the distribution 
of the squeezing between conjugate variables. The squeezed vacuum state occurs when 
a = 0: 



The limit of high squeezing occurs when r >> 1 , where the state \z > is highly localised in 
momentum space. 

Single-mode squeezed operators do not conserve momentum, since they describe the 
creation of particle pairs with momentum k . Two-mode squeezed operators, however, 
describe the creation and annihilation of two particles (waves) with equal and opposite 
momenta. A two-mode squeeze operator has the form [32] 



where a± and a± are annihilation and creation operators for the two modes, respectively. 

Consider two conjugate operators p and q , with variances Ap and Aq . In the 
squeezed-state formalism, one may construct states such that Ap and Aq are equal, taking 
the minimum value possible. The name 'squeezed' refers to the fact that the variance of one 
variable in a conjugate pair can go below the minimum allowed by the uncertainty principle 
(the squeezed variable), while the variance of the conjugate variable can exceed the min- 
imum value allowed (the superfluctuant variable) [25,33,34]. The superfluctuant variable 



z> = \0,z> = S(z)\0> 



(5.7) 




(5.8) 
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is amplified by the squeezing process, and so becomes possible to observe macroscopically, 
while the subfluctuant variable is squeezed and becomes unobservable. In particle pro- 
duction, whether by black holes or in cosmology, the number operator is a superfluctuant 
variable, while the phase is squeezed. 

6. Analytic continuation and the large-squeezing limit 

We shall see here for the black- hole evaporation problem that, when one rotates the 
time-separation T at infinity: T — > |T| exp(— id) into the complex by a very small angle 
8 > , one arrives at a very highly-squeezed quantum state. There is no information-loss 
paradox associated with the relic Hawking radiation, as such a state is a pure state. It 
is also important to state that we do not take the \T\ — > oo limit. However, one must 
understand that the observation time at infinity by far exceeds the dynamical collapse 
time-scale, which is of order irMj [23]. We now repeat Eq.(2.19): 



{A,}; t\ = N e~ lMlT ]| tf(A,-;T) 



(6.1) 



and then define 



{A,}; T] = N e~ lMlT Y[ 2i sin(^T) #(A,-;T) 



= Ne 



-iMi 



(Akj) kj |A/ cot(^T) 



IlT J] exp 

3 

= Ne X p(iS$L B [{A j };T]). 

We further define the functions <j>j ( |T| , 8 ) and Tj ( \T\ , 8 ) by 

<f)j(\T\,6) = — kj\T\ cos 8 , 
tanhrj( \T\ ,8) = exp(— 2kj \T\ sin 8 ) , 

whence 

exp(— 2rj) = tanh(/cj |T| sin 8 ) 
From Eqs. (6.3-5), one can rewrite Eq.(6.2) in the form 



{Aj};\T\,8 



= Ne 



-iMj|T|cos5 C -Mj |T| sin 5 



n 



exp 



■-(Afe^-f — 



1 + e 2i ^ tanh: 



e 2 *^ tanhr. 



(6.2) 



(6.3) 
(6.4) 

(6.5) 



(6.6) 



On comparing with Sec. 5, we recognise Eq.(6.6) as the coordinate-space representation of 
a quantum-mechanical squeezed state [35,36], with rj( \T\ ,8) the squeeze parameter and 
4>j( \T\,8) the squeeze angle. The evolution of the squeezed state is taken into account by 
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the |T|-dependence in rj and in cj)j , which are in general both complicated functions of 
time. 

We now define 

ej = kj \T\ sin8 , |T|) = i + sm2 ( fc J l T D _ ( 6 7 ) 



sinh €j 



Then 

|$[{^-};|T|,5]| 2 = |AT| 2 e - 2M ^ T l sin5 J] exp 



3 



/(*j,£j.m) 



, (6.8) 



and, from Eqs. (6.5,7): 



~ e" 2 ^ , ej<l , (6.9) 



corresponding to ^> 1 , which is the limit of high squeezing. We discuss the form of the 
normalisation in another paper [9]. 

Eq.(6.8) describes a Gaussian non-stationary process in which the variance is an oscil- 
latory function of time. Rather than dealing with travelling waves, one now has standing 
bosonic waves, where the amplitudes for left- and right-moving waves are large and almost 
equal - this is similar to the inflationary-cosmology scenario [17]. One consequence of the 
high-squeezing behaviour is that the variance for the amplitudes {xj} is large, so that 
there are large statistical deviations of the observable power spectrum from its expected 
value. This is just a manifestation of the Uncertainty Principle. 

In the squeezed- state formalism, the high-squeezing limit rj ^> 1 may be regarded 
as the classical limit. For example, in this sense, in the case of black-hole evaporation, 
the final state of the remnant particle flux becomes more classical (more WKB) in the 
limit 5^0. In this limit, one can effectively consider the final perturbations as being 
represented by a classical probability distribution [17,33,37]. As in the inflationary sce- 
nario in cosmology, the perturbations on the spherically-symmetric black-hole background 
space-time, of quantum-mechanical origin, cannot be distinguished from classical stochas- 
tic perturbations, without the need for an environment for decoherence. There is also a 
correspondence between the initial conditions for the perturbations in the black-hole and 
in the cosmological cases. In cosmology, the assumption is that, at some early 'time' just 
prior to inflation, the modes are in their adiabatic ground state. A similar qualitative 
statement can be made in the black-hole example, provided that the pre-collapse initial 
data were diffuse, slowly-moving and spherically symmetric. 

One further consequence follows, provided that ej is small (as above). Then, one finds 
for the probability distribution Eq.(6.8) that, as 6 — > + , 

oo 

|$[{^-};|T|,(J]| 2 ~ |iV| 2 J] J] exp[-(A^ n )c^|A sn , mP | 2 ] , (6.10) 

slmP n=l 
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where we have used the approximation sinh ej ~ ej for small 6j , and the identities 

1 

S (x) = — lim 



and 

*[/(*)] 



7T e^O (e 2 + X 2 ) 

E5(x — Xi) 



where Xi are zeros of f(x) and u n = mr/\T\ , Au n = (iv n+ i — u> n ) . We have also used 
the fact that kj — > and that kj \Aj\ — > as kj — > . In practice, the product over n 
should be cut off at some large n max , such that uj nmax = Mj . 

Further investigation of the derivation of Eq.(6.10) indicates that, in the limit of high 
squeezing, the random variable <f>j associated with the final state is squeezed to discrete 
values, independently of the quantum numbers {simP} [9]. Note that it is only the 
squeeze phases {4>j} of the (standing- wave) perturbations which are fixed and correlated 
in the high-squeezing limit. 

For comparison, in inflationary cosmology, the oscillation phases of standing waves 
have fixed values, giving rise to zeros in the power spectrum, which are characteristic of 
the CMBR. The power spectrum of cosmological perturbations in the present universe is 
not a smooth function of frequency. The standing- wave pattern, due to squeezing, induces 
oscillations in the power spectrum. This in turn produces Sakharov oscillations [37,38], 
due to metric and scalar perturbations in the distribution of higher-order multipoles of 
the angular correlation function for the temperature anisotropies [21,39] in the CMBR, 
for all perturbations at a given time whose wavelength is comparable with or greater than 
the Hubble radius defined for that time. That is, the peaks and troughs of the angular 
power spectrum have a close relationship with the maxima and minima of the metric power 
spectrum. For long wavelengths, the power spectrum does become smoother. 

7. Conclusion 

In this Letter, we have illustrated many aspects of the quantum boundary- value for- 
mulation, for linearised bosonic fields (spins s = , 1 , 2) propagating in the space-time of 
an evaporating black hole. When the Lorentzian proper-time separation T between the 
initial and final space-like hypersurfaces, as measured at spatial infinity, is deformed into 
the lower complex T-plane, and when the perturbations are initially weak, one obtains a 
quantum-mechanical squeezed- state formalism. The large-squeezing limit is equivalent to 
the WKB limit, corresponding to an infinitesimal angle 5 of rotation of T into the 
lower-half complex plane. 

Since the final squeezed state is a pure state, there is no information-loss paradox 
as a result of the Feynman + ie prescription we have adopted. Our complex approach is 
new and differs from Grischchuk's original application of squeezed states to black holes. 
However, as in the cosmological scenario, so the bosonic perturbations on the black-hole 
background can be regarded as a stochastic collection of standing waves, rather than as 
travelling waves, in the high-squeezing limit. This leads to the prediction of peaks in the 
power spectrum of the relic black-hole radiation, analogous to the Sakharov oscillations in 
the CMBR. 
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